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Abstract. Broad absorption lines (BALs) in quasar spectra identify high velocity 
outflows that likely exist in all quasars and could play a major role in feedback to 
galaxy evolution. Studying the variability in these BALs can help us understand the 
structure, evolution, and basic physical properties of these outflows. We are conducting 
a BAL monitoring program, which so far includes 163 spectra of 24 luminous quasars, 
covering time-scales from ~1 week to 8 years in the quasar rest-frame. We investigate 
changes in both the CIV /I1550 and SilV /11400 BALs, and we report here on some of 
the results from this program. 



1. Introduction 

High-velocity wind-like outflows are an important part of the quasar system and a po- 
tential contributor to feedback to the host galaxy. Many properties of quasar outflows 
are still poorly understood, including their location and three-dimensional structure. 
One way to obtain constraints on quasar outflows is to study the variability in their ab- 
sorption lines. Variability on shorter time-scales can place constraints on the distance 
of the absorbing material from the central supermassive black hole (SMBH). Measure- 
ments of variability on longer (multi-year) time-scales provide insight into the homo- 
geneity and stability of the outflowing gas. Overall, results of variability studies provide 
information on the size, kinematics, and internal makeup of sub-structures within the 
outflows. 

Broad absorption lines (BALs) are the most prominent signatures of accretion disk 
outflows seen in quasar spectra. We are conducting a BAL monitorin g program for a 



sample of 24 luminous BAL quasars at z= 1.2-2.9 from iBarlowl ( 19931) . We have been 



re-observing these quasars at the MDM Observatory 2.4-m Hiltner and KPNO 2. 1-m 
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telescopes, and we also include spectra from SDSS when available. We currently have 
in total 163 spectra of these 24 quasars. 

To identify BAL variability, we first define the regions of BAL absorption, based 
on the definition of the 'balnicity index' (BI), i.e. they must contain cont iguous absorp- 
tion t hat reaches > 10 per cent below the continuum across >2000 km/s (IWeymann et al. 
Il99lh . We then used visual inspection to identify velocity intervals with a width of at 
least 1200 km/s that varied. Using the average flux and associated error within each 
candidate variability interval, we place an error on the flux differences between the two 
spectra. We include all intervals of variability where the flux differences are at least 4cr. 



2. CIV and SilV BAL Variability 

2.1. Trends in the CIV Short-Term and Long-Term Data 



The first results of our BAL variability study were reported in lCapellupo et all (1201 la : 



hereafter, Paper 1), where we focused on variability in just the CIV /U550 BALs in 
two different time intervals: 4-9 months (short-term) and 3.8-7.7 yr (long-term) in the 
quasar rest-frame. We found both a higher incidence of variability (65%, or 15/23, of 
the quasars versus 39%, or 7/18) and a larger typical change in strength in the long-term 
data. Th e variability occurs typically in only portions of the BAL troughs (e.g. Fig. [3j 
see also, iGibson et al.ll2010l ). In Fig. [TJ we plot the incidence of CIV BAL absorption 
and variability versus velocity. The BAL components at higher outflow velocities are 
more likely to vary than those at lower velocities. We also find in Paper 1 that weaker 
BALs are more likely to vary than stronger BALs. 

2.2. Multi-Epoch Monitoring of SilV and CIV Variability 

In ICapellupo et all d2011al : hereafter, Paper 2), we directly compare the variability of 



SilV /11400 BALs to CIV in the long-term dataset from Paper 1. C and Si have different 
abundanc es, if solar abu ndances are assumed, and they have different ionization proper- 



ties (e.g. lHamann et al.ll2008l) . By examining if CIV and SilV have different variability 



properties, and how they differ, coupled with these differences in abundances and ion- 
ization properties, we can gain new insight into the cause(s) of BAL variability. 

SilV BALs are more likely to vary than CIV BALs. For example, when looking 
at flow speeds >-20 000 km/s, 47 per cent of the quasars in our sample exhibited 
SilV variability while 31 per cent exhibited CIV variability. Variability in SilV can 
occur without corresponding changes in CIV at the same velocities. For ~50 per cent 
of the variable SilV regions, there was no corresponding CIV variability at the same 
velocities. However, we have only one tentative case where changes in CIV are not 
matched by SilV. At BAL velocities where both CIV and SilV varied, the changes 
always occurred in the same sense (with both lines either getting stronger or weaker). 

We now have up to 13 epochs of data per object for our sample of 24 BAL quasars. 
With all the observing epochs included, the fraction of quasars with CIV BAL variabil- 
ity at any velocity increases from the value of 65% found in Paper 1 to 88% for the 
same sample of quasars. This increase was caused by variations missed in the 2-epoch 
comparisons in Paper 1. We also find that BAL changes at different velocities in the 
same ion almost always occurred in the same sense. We found just 3 cases that show 
evidence for one CIV BAL weakening while another strengthens within the same ob- 
ject. The multi-epoch data also show that the BAL changes across ~ 1 week to 8 years 
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Figure 1 . The top two panels show the number of occurrences of CIV BAL ab- 
sorption and BAL variability versus velocity. The third panel is the second panel 
divided by the first, with l<x error bars overplotted. 



in the rest-frame were not generally mono tonic (see also, iGibson et al.ll20ldh . Thus, 
the characteristic time-scale for significant line variations, and (perhaps) for structural 
changes in the outflows, is less than a few years. Furthermore, with more epochs added, 
we still do not find clear evidence for acceleration or deceleration in the BAL outflows. 



2.3. Time-Scales of Variability 

We take a closer look at the time-scales of variability in Capellup o et al.1 ( 20121 ; here- 
after, Paper 3). In Spring 2010, we re-observed a subsample of our BAL quasars mul- 
tiple times over rest-frame time intervals of At ~ 1 week to 1 month to augment our 
temporal sampling at these short time-scales. Variability results on the shortest time- 
scales are important for putting constraints on the location and sizes of the outflows. 

The multi-epoch nature of our dataset, spanning time intervals from 0.02 to 8.7 yr 
in the rest frame, allows us to investigate the probability of BAL variability versus At. 
We examined CIV BAL absorption at all measured velocities in our entire data sample. 
We compared each pair of observations for each quasar and counted the occurrences 
of CIV BAL variability, using our definition of BAL variability defined in Paper 1 (see 
Section 1 above). We then calculated a probability by dividing the number of occur- 
rences of variability by the number of measurements, where a pair of observations is 
one measurement, in logarithmic bins of At. We plot this measured probability of de- 
tecting CIV BAL variability versus At in Fig. |2] The middle and bottom panels display 
the number of measurements and the number of quasars, respectively, that contribute to 
each At bin. Each quasar can contribute multiple times to each bin. 
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Figure 2. The top panel displays the fraction of measurements where there was 
CIV BAL variability versus the time interval between the two observations. The 
middle and bottom panels show the number of measurements and the number of 
quasars, respectively, that contribute to each bin. 



In our dataset, we found variability down to a At of 8 days in 1246-0542 between 
the 2010.13 and 2010.20 observations. Fig. |3] shows the CIV region of the spectra 
for these two observations as bold and thin lines, respectively. As indicated by shaded 
bars in Fig. |3j there are two separate velocity intervals of variability between the two 
spectra, and the flux difference is at least 8cr in each interval. We also overplot an earlier 
(2007.04; thin dashed curve) and a later (2010.35; bold dashed curve) epoch in Fig. |3] 
which shows that there was variability in these intervals in other epochs as well. 



3. Conclusions 

Coordinated variabilities between absorption regions at different velocities in individ- 
ual quasars seems to favor changing ionization of the outflowing gas as the cause of the 
observed BAL variability. Many of our results are consistent with this scenario, includ- 
ing our findings that SilV, which has a lower optical depth and is therefore less likely to 
be saturated than CIV, is more likely to vary than CIV. However, variability in limited 
portions of broad troughs fits naturally in a scenario where movements of individual 
clouds, or substructures in the flow, across our lines-of-sight cause the absorption to 
vary. This scenario is also consistent with the main results of Paper 2 (Section 2.2 
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Figure 3. Smoothed MDM spectra of 1246-0542, showing the CIV absorption 
region. The bold and thin solid curves show the two epochs separated by 8 days 
in the rest-frame, and the shaded regions indicate the velocity intervals that varied 
between these two epochs. The formal lcr errors are shown near the bottom. 

above) if SilV has a smaller covering fraction than CIV. The actual situation may be a 
complex mixture of changing ionization and cloud movements (Paper 2). 

If the short time-scale variations in 1246-0542 were due to cloud movements 
and a ~10% change in line strength over 8 days indicates the cloud crossed 10% of 
the continuum source in that time, then the transverse speed is ~-27 000 km/s. If 
we assume the transverse speed is related to the Keplarian rotation speed, then the 
outflowing gas is located <scl pc from the central SMBH, probably inside the radius of 
the broad emission-line region (Paper 3). 
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